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摘 要 : 反 演 成 像 是 综合 孔径 微波 辐射 计 的 一 项 关键 内 容 , 但 反 演 过 程 是 欠 定 的 ,不 能 提供 一 个 确定 解 . 针对 综合 
孔径 辐射 计 中 常用 反 演 算法 :G 矩阵 模型 法 存在 较 大 反 演 误差 的 问题 ,提出 了 一 种 改进 的 G 矩阵 模型 法 应 用 于 
FPIR 系统 中 . 仿真 结果 表明 ,与 传统 的 G 和 矩阵 模型 法 相 比 ,改进 的 G 和 矩阵 模型 法 能 有 效 地 降低 FPIR 系统 图 像 反 
演 误差 ,以 获取 高 精确 度 的 观测 场景 的 亮 温 分 布 满足 FPIR 系统 探测 海面 风 场 .土壤 湿度 等 应 用 需求 . 
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Application of Improved G Matrix Model in Full 
Polarization Interferometric Radiometer 
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Abstract: The inverse imaging is a key issue for microwave interferometric radiometers, but the 
inverse process is iposed and doesn't provide a unique and stable solution. For the problem that 
there exists a large reconstruction bias with G matrix model used generally in the inversion 
algorithm, an improved G matrix model applied in FPIR system was presented. The simulation 
results show that, compared with conventional G matrix model, the improved G matrix model can 
greatly reduce the reconstruction bias in the FPIR system. The improved G matrix model can 
obtain more accurately the brightness temperature distribution of the observed scene, so as to 
meet requirements of measuring wind field of the ocean surface and soil moisture for FPIR. 
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综合 孔径 微波 辐射 计 是 从 20 世纪 80 年 代 后 期 soil moisture and ocean salinity) 卫 星 的 发 射 成 功 为 
开始 兴起 的 一 项 被 动 微 波 遥 感 新 技术 . 伴随 着 对 干 标志 ,综合 孔径 微波 辐射 计 已 经 逐渐 成 熟 并 走向 了 
涉 式 综合 孔径 成 像 理 论 研究 的 深入 以 及 数字 技术 的 实际 应 用 5. 一 维 综合 孔径 辐射 计 具 有 体积 小 、 结 
快速 发 展 , 以 ESTARD (electronically scanned 为 简单 、 功 耗 重 量 低 等 优点 ,但 是 传统 的 一 维 综合 
thinned array radiometer) 的 试飞 和 SMOS (the 径 辐 射 计 只 能 获取 单 极 化 或 者 双 极 化 的 信息 ,无 法 
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获取 全 极 化 信息 ,因此 一 些 重 要 的 地 球 物理 参数 ,如 
海洋 表面 的 风 场 信息 等 ,就 无 法 通过 传统 的 一 维 综 
合 孔 径 辐 射 计 获 得 . 全 极 化 参量 微波 辐射 计 是 另 一 
项 被 动 微波 遥感 的 新 技术 ,能 够 获得 全 极 化 信息 , 主 
要 应 用 于 探测 海面 的 风 场 信息 . 

一 维 全 极 化 综合 孔径 辐射 计 (full polarization 
interferometric radiometer，FPIR) 是 一 套 高 分 辩 
率 、 轻 量化 、. 低 功 耗 的 X 波 段 微 波 辐射 计 " 一 . 它 将 
干涉 式 综 合 孔径 成 像 技术 和 全 极 化 参量 信息 获取 结 
合 , 与 ESTAR 等 传统 的 一 维 综合 孔径 辐射 计 不 同 
之 处 是 全 极 化 信息 获取 ,圆锥 面 天 线 波 束 和 两 点 定 
标 等 . FPIR 的 应 用 目标 之 一 是 搭载 在 欧 空 局 
SMOS 业务 卫星 上 ,与 MIRASCmicrowave imaging 
radiometer with aperture synthesis) 一 起 工作 ,形成 
一 套 双 频 综 合 孔径 微波 辐射 计 系 统 . FPIR 产生 X 
波段 的 全 极 化 亮 温 图 像 ,为 MIRAS 提供 海面 风 场 
和 陆地 植被 的 信息 ,从 而 提高 了 MIRAS 海水 盐 度 
和 土壤 湿度 的 反 演 精度 . 

反 演 成 像 是 综合 孔径 辐射 计 的 一 项 关键 内 容 . 
在 综合 孔径 辐射 计 中 ,常用 的 反 演算 法 是 侍 里 叶 变 换 
法 和 G 和 矩阵 模型 法 ,其 中 G 和 矩阵 模型 法 将 亮 温 图 
像 反 演 与 误差 定 标 相 结合 ,在 综合 孔径 辐射 计 成 像 中 
得 到 广泛 应 用 . 但 是 即使 在 没有 模型 误差 和 辐射 噪 
声 情况 下 ,G 矩阵 模型 法 依然 存在 着 较 大 的 反 演 误 
差 . 为 了 降低 图 像 反 演 误差 , 本文 提出 了 一 种 改进 的 
G 矩阵 模型 法 ,并 基于 FPIR 系统 ,进行 了 仿真 . 


1 FPIR 系统 的 成 像 原 理 


FPIR 的 两 单元 原理 框图 如 图 1 所 示 . 在 顺 轨 
方向 ,天 线 方向 图 为 真实 孔径 线 阵 形成 的 窄 波 束 ;在 
交 轨 方向 ,天 线 方向 图 由 通常 的 局 形变 成 圆锥 面 ,每 
个 像素 的 入 射 角 a 相同 . 

复 可 视 度 画 数 V(w; ) 与 辐射 亮 温 T,(6) 之 间 存 


坐标 系 天 线 i 天 线 平面 (z=0) 
- 接收 机 上 [ADCH -=| 数 


强 水 昔 上 
a 


| [ 柜 & 机 一 [ADd] 


EF 面 (z= 有 ) 


天 线 足迹 ~ 
y 
图 ] FPIR 系统 的 两 单元 原理 框图 


Fig.1 The schematic block diagram of a two-element FPIR system 


在 如 下 关系 


”sin 0 个 (E) TT 
Va) =| Ep, Ft, Xx 
| 一 sin 0 sin20 一 e 
EC Ye)exp(—j2rust) dé (1) 
区 


式 中 局 二 dj; /和 为 天 线 i 和 j 之 间 归 一 化 的 基线 长 
度 ;E 二 sin 0sin $ 为 方向 余弦 ,0 和 y 为 球面 坐标 ; fo 
为 中 心 频率 ; T, 为 接收 机 的 物理 温度 ; F (6) 和 
FX (各 分别 为 天 线 单元 i 和 j 所 对 应 的 归 一 化 的 电 
压 方向 图 (* 表示 共 斩 );r， (Ar) 称 为 条 纹 洗涤 
(fringe washing) 落 数 ,在 信道 带宽 足够 小 的 情况 
下 ,有 rxr(Ar) 守 1. 为 了 简化 , 今 5 一 6/sin 0 一 sin 多， 
则 式 (1) 变 为 


i , 
Vlws) = | 2 Fi(O F(X 
-1 Mie : ? 
RE 人 exp(— ji2rust) ds. ( 2) 


将 式 (2) 与 传统 的 一 维 综合 孔径 辐射 可 视 度 画 数 与 
辐射 亮 温 之 间 的 公式 比较 ,形式 上 的 不 同 之 处 是 变 
量 £ 二 fsin 9, 但 实际 上 地 球 物理 参数 的 反 演 却 得 到 
简化 . 
将 式 (2) 进 行 离散 化 得 
Vontvxi = CoNHDxxM wxl， (3) 
天 自 | 吉 :日 td s 
其 中 修正 的 辐射 亮 温 工 (9) 一 被 离散 为 
M 项 ,可 视 度 画 数 V(w; ) 被 离散 化 为 2N 十 ! 项. 通 
常情 况 下 选取 的 亮 温 反 演 点 数 M 均 大 于 可 视 度 画 
数 采样 点 数 2N 十 1 的 3 倍 , 所 以 式 (3) 是 欠 定 . 假 
设 辐射 亮 温 的 工作 空间 为 如, 可 视 度 函 数 的 工作 区 
间 为 C, 则 根据 最 小 二 乘 准则 


min|v 一 GT|:. (4) 
TEB 
求解 其 最 小 范 数 解 ,得 
T= GtVy = G(GG") Vy. (5) 


式 中 :个 为 反 演 的 亮 温 ;G1 = 二 GH(GG?) 1! 为 广义 北 
和 矩阵 . 

为 了 降低 空间 频率 上 有 限 的 采样 覆盖 所 导致 的 
吉 布 斯 振荡 现象 ,对 可 视 度 采样 进行 加 窗 处 理 , 即 

个 . = UWUT. (6 ) 

式 中 :U 为 离散 伟 里 叶 变 换 ;W 为 窗 画 数 . 常用 的 窗 
画 数 有 : 汉 明 窗 、 汉 宁 窗 、 布 莱克 曼 窗 、 契 比 雪夫 窗 
等 . 值得 注意 的 是 ,在 分 析 反 演 图 像 与 真实 图 像 的 
误差 时 ,需要 将 加 窗 后 反 演 结果 个 , 与 真实 图 像 加 窗 
后 结果 T. 王 U8WUT 进行 比较 . 
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2 改进 的 G 和 矩阵 模型 法 


利用 G 和 矩阵 模型 法 重 构 观 测 目标 亮 温 时 发 现 ， 
即使 在 不 存在 G 和 矩阵 模型 误差 和 辐射 噪声 情况 下 ， 
依然 存在 图 像 反 演 误差 . 而 图 像 反 演 误 差 只 能 在 图 
像 反 演 过 程 中 去 除 或 者 降低 ,而 不 能 利用 设备 定 标 
方法 进行 去 除 或 者 降低 . E. Anterrieu 等 中 在 分 析 
MIRAS 系统 反 演 过 程 中 发 现 ,图 像 反 演 误差 不 仅 
来 源 于 观测 目标 的 高 频 分 量 的 缺失 ,还 与 参考 辐射 
计 天 线 方向 图 和 目标 无 温 荆 视 场 外 亮 温 对 视 场 内 亮 
温 反 演 的 影响 有 关 . 而 对 于 FPIR 系统 ,研究 表明 : 
图 像 反 演 误差 不 仅 与 无 混 荆 视 场 外 亮 温 对 视 场 内 亮 
温 反 演 的 影响 有 关 , 而 且 与 单元 天 线 方 向 图 误差 有 
关 ,特别 是 单元 天 线 方向 图 幅度 的 起 伏 及 天 线 方向 
图 之 间 的 不 一 致 . 

为 了 降低 FPIR 系统 图 像 反 演 误差 ,提出 了 一 
种 新 的 降低 系统 反 演 误差 的 方法 , 即 :改进 的 G 算 
阵 模型 法 . 这 种 方法 是 基于 E. Anterrieu 等 号 提出 
的 算法 ,引入 一 个 尽 可 能 接近 观测 场景 真实 亮 温 了 
的 亮 温 图 像 工 , 则 


8Vy 一 站 一 站 一 了 一 GT. (7) 
则 反 演 亮 温 变 为 
个 二 G+ (YY 一 GT) 十 工 (8) 


为 了 得 到 亮 温 图 像 了 ,最 简单 的 方法 就 是 根据 
海岸 线 等 先 验 信息 ,将 观测 场景 划分 为 陆地 区 域 和 
海洋 区 域 . 假定 i 是 陆地 区 域 的 亮 温 , 设 定 为 恒定 
的 ;to 是 海洋 区 域 的 亮 温 , 设 定 为 恒定 的 . 亮 温 图 像 
了 由 ti 和 zio 构成 ,T 可 以 表示 为 

T= Dota (9) 
式 中 DD 由 Di 和 Do 组 成 ,Di 和 Do 分 别 为 与 陆地 
和 海洋 区 域 对 应 的 二 元 掩 膜 (binary mask); 


1 一 [万 i6|: 
根据 最 小 二 乘 准则 ,得 
min|y 一 VCz) |:. (10) 
1{ER” 
求解 可 得 
1 一 [(GD)aCGD)] (GD)'V. (11) 
以 上 方法 存在 一 个 明显 的 缺点 :存在 两 个 不 同 


的 最 小 二 乘 准则 , 即 : 式 (4) 和 式 (10). 由 于 两 个 最 
小 二 乘 准则 不 能 同时 满足 ,所 以 为 了 解决 这 个 问题 ， 
提出 了 一 种 新 的 算法 计算 人 工 图 像 工 这 种 算法 是 
基于 Ali Khazaal 提出 的 算法 ,利用 最 小 二 乘 准 
则 , 即 


min min|V VCO) — GIT— TO0)]|:. (12) 
2 TE 


iER” 
则 反 演 亮 温 为 
多 一 GT (YY 一 GT) 十 工 (13) 
其 中 /上 满足 
min|V— VO) 


LERP 


min|lV— VG) 一 GG+ [Vo VG): = 


LE RP 


min| (I— GG EV — VO) |:. (14) 


求解 可 得 

1 一 [COPGCD)aCPGD)] 一 (CPGD)aPV， (15) 
其 中 P==I 一 GG+. 当 P= 工 时 式 (15) 与 式 (11) 相 
同 . 为 了 区 别 ,将 式 (11) 方 法 命名 为 改进 的 G 矩阵 
方法 1, 式 (15) 方 法 命名 为 改进 的 G 和 矩阵 方法 2. 

上 述 的 改进 G 和 矩阵 模型 法 ,是 基于 E，、Anter- 
rieu 和 Ali Khazaal 提出 的 算法 . 但 与 E. Anterrieu 
和 Ali Khazaal 提出 的 算法 相 比 ,应 用 于 FPIR 的 方 
法 有 以 下 不 同 之 处 . 首先 ,与 MIRAS 系统 相 比 ， 
FPIR 系统 反 演 过 程 简 单 得 多 ,所 以 计算 效率 更 高 ， 
所 费 存 储量 更 小 . 其 次 ,FPIR 形成 圆锥 面 天 线 波 
束 , 在 地 物 平 面 的 每 个 像素 的 入 射 角 相 同 ,而 MIR- 
AS 中 在 地 物 平 面 的 不 同 像素 的 入 射 角 不 同 , 所 以 
FPIR 中 对 工 建 模 前 的 先 验 知识 不 需要 包括 入 射 角 
信息 . 同时 ,FPIR 中 工 仅 由 和 和 to 构成 ,不 受 天 空 
等 其 他 参数 的 影响 ,因此 建 模 更 加 简单 . 最 后 ,应 用 
于 FPIR 的 这 种 方法 是 与 传统 G 和 矩阵 模型 法 相 结 
合 , 有 效 地 去 除 传统 G 矩阵 模型 法 中 的 反 演 误差 ， 
而 E. Anterrieu 和 Ali Khazaal 提出 的 方法 是 与 解 
析 和 矩阵 方法 相合 . 


3 仿 真 


为 了 验证 上 述 反 演 算法 的 正确 性 ,进行 了 仿真 
分 析 . 选择 仿真 对 象 为 FPIR 系统 工作 在 H 极 化 状 
态 ,天 线 单元 数量 /一 9, 最 短 基线 为 Ad 一 0 635N， 
最 长 基线 为 25Ad ,忽略 空间 去 相关 效应 的 影响 综 
合 孔 径 方 向 上 单元 天 线 功 率 方向 图 如 图 2 所 示 ,为 
章 室 中 真实 测量 结果 . 

仿真 原始 图 像 如 图 3 所 示 . 图 3 为 天 线 平面 内 
亮 温 ,来 源 于 我 国 渤海 地 区 的 光学 卫星 遥感 图 像 . 

3 中 ,垂直 方向 为 真实 孔径 方向 ,水 平方 向 为 
综合 孔径 方向 ， 对 于 一 维 综合 孔径 辐射 计 , 单 次 测 
可 获得 视 场 内 的 一 个 条 带 的 图 像 ,通过 平台 运动 
行 推 扫 从 而 获得 二 维 图 像 . 


GIT— TQ)]|: = 


内 


洋 甩 
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示 , 所 使 用 的 窗 画 数 为 汉 宁 窗 . 为 了 便于 比较 ,无 混 
鸽 视 场 内 原始 亮 温 图 像 加 窗 后 结果 如 图 4(a) 所 示 . 
所 240 
其 230 
国 220 
Ey 210 . 
200 
1 1 1 | 证 190 到 
-60 -40 -20 0 20 40 60 到 180 
方位 角 /。 ) a |: 
图 2 综合 孔径 方向 天 线 单元 方向 图 8 
Fig.2 Synthetic aperture patterns of 9 antennas 6 


分 别 利 用 传统 的 G 矩阵 法 改进 的 G 和 矩阵 方法 
1 改进 的 G 矩阵 方法 2 得 到 了 其 对 应 的 无 混 亚 视 
场 内 加 窗 后 反 演 的 结果 ,如 图 4(b)、4(Cc) 和 4(d) 所 
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Fig.3 The actual image 
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Fig.4 The reconstructed images (brightness temperature) 


为 了 定量 地 分 析 加 窗 后 成 像 结 果 与 原始 图 像 的 
误差 AT, 二 T, 一 必 , ,引入 偏差 和 标准 差 . 偏差 和 标 
准 差 的 计算 公式 分 别 为 


= >) (个 一 


AFOV 


= TT 
式 中 :AT, 为 图 像 中 某 行 即 某 条 带 的 偏差 ;5, 为 某 
行 的 标准 差 ; zsrov 为 某 行 无 混 赤 区 域内 的 像素 总 
数 ; 们 为 某 行 第 & 个 像素 的 反 演 亮 温 ; T, 为 与 个 ， 
对 应 的 原始 亮 温 . 

对 于 图 4 的 反 演 结果 ,偏差 和 标准 差 对 比 结 


Ti; ) /narov. (16) 


— AT,)’/nagoy ， GCL? 


如 图 5 和 图 6 所 示 . 图 中 横 轴 为 FPIR 顺 轨 推 扫 时 ， 
按时 间 顺 序 排列 的 观测 条 带 编 号 , 纵 轴 分 别 为 不 同 
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真实 孔径 方向 推 扫 的 条 带 /个 
图 5 偏差 对 比 


Fig.5 The bias contrast 
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有 反 演算 法 得 到 的 反 演 结果 的 偏差 和 标准 差 . 从 图 中 
可 以 看 出 ,与 传统 的 G 和 矩阵 模型 法 相 比 ,改进 的 G 
答 阵 方法 1 能 够 有 效 地 降低 残留 的 图 像 反 演 误差 ， 
但 在 某 些 区 域 反 演 结果 会 出 现 恶 化 ,这 是 由 于 两 个 
最 小 平方 准则 不 能 同时 满足 导致 的 . 而 改进 的 G 和 矩 
阵 方法 2 克服 了 这 个 缺陷 ,有 效 地 降低 了 残留 的 图 
像 反 演 误差 . 值得 注意 的 是 , 反 演 误差 与 地 物 目标 
的 特征 相关 , 随 着 观测 目标 的 变化 , 反 演 误差 也 产生 
相应 的 变化 . 
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真实 孔径 方向 推 扫 的 条 带 /个 
图 6 标准 差 对 比 


Fig.6 The standard deviation contrast 


4 结 论 


FPIR 是 一 套 高 分 辨 率 、 轻 量化 、 低 功 耗 的 X 波 
段 一 维 综合 孔径 微波 辐射 计 . 它 与 传统 一 维 综合 孔 
径 辐 射 计 的 不 同 之 处 在 于 全 极 化 信息 获取 ,圆锥 面 
天 线 波束 和 两 点 定 标 等 . 反 演 成 像 是 干涉 式 综合 孔 
径 辐 射 计 的 一 项 关键 内 容 . 因为 辐射 计 干 涉 输出 的 
结果 是 可 视 度 本 数 采样 ,所 以 如 何 准 确 有 效 地 从 可 
视 度 画 数 的 采样 数据 转化 到 其 空间 域 的 辐射 图 像 ， 
是 成 像 反 演算 法 所 要 解决 的 问题 . 但 由 于 反 演 过 程 
是 欠 定 的 ,无 法 得 到 一 个 准确 解 ,即使 在 没有 模型 误 
差 和 辐射 噪声 情况 下 ,依然 存在 图 像 反 演 误差 . 

为 了 降低 FPIR 系统 图 像 反 演 误 差 ,提出 了 一 
种 改进 的 G 和 矩阵 模型 法 . 仿真 结果 表明 ,与 传统 的 
G 和 矩阵 模型 法 相 比 ,改进 的 G 矩阵 模型 法 能 有 效 地 
降低 图 像 反 演 误差 ,证 明了 改进 的 G 矩阵 模型 法 在 
FPIR 中 应 用 的 有 效 性 . 同时 ,改进 的 G 矩阵 模型 法 
也 可 以 应 用 于 其 他 的 一 维 综合 孔径 微波 辐射 计 中 . 
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